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Abstract 

We report the amplitude scaling behavior of Frenkel exciton chains with nearest- 
neighbor correlated off-diagonal random interactions. The band center spectrum 
and its localization properties are investigated through the integrated density of 
states and the inverse localization length. The correlated random interactions are 
produced through a binary sequence similar to the interactions in spin glass chains. 
We produced sets of data with different interaction strength and "wrong" sign 
concentrations that collapsed after scaling to the predictions of a theory developed 
earlier for Dirac fermions with random- varying mass. We found good agreement 
as the energy approaches the band center for a wide range of concentrations. We 
have also established the concentration dependence of the lowest order expansion 
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coefficient of the scaling amplitudes for the correlated case. The correlation causes 
unusual behavior of the spectra, i.e., deviations from the Dyson-type singularity. 



1 Introduction 

Low- dimensional disordered excitonic systems have generated renewed interest re- 
cently. For example, the possibility of using optically active polymers in various pho- 
tonic applications such as flat panel displays and light emitting diodes has led to 
increased interest in their electronic and optical properties [TH21H] . We study here nu- 
merically the dynamics of Frenkel excitons near the band center in a system with 
nearest-neighbor off-diagonal random interactions. This study is, we believe, the first 
to treat correlated disorder in the off-diagonal dipolar interactions in an exciton sys- 
tem which focuses on the behavior near band center E = 0. However, an analysis of a 
similar problem in the correlated electronic system has appeared [Sj but the authors 
considered correlations those were different from ours, and the interactions were not of 
the dipolar type [3] . Many efforts have also been spent on the same off-diagonal tight- 
binding system using supersymmetric methods (SUSY) [B] where the interactions are 
formally similar to our dipolar interactions but the authors eventually worked with 
continuous variables whereas ours are discrete. Later, correlated disorder, in partic- 
ular the exponential type, was included for the same tight binding model using the 
SUSY [7118] . Our results can thus be compared with what has been developed for the 
SUSY methods. This study is important because most realistic random disordered 
systems have nonlocal correlations. 

The nature of the excitations at the band center of random off-diagonal exciton models 
has been discussed for a long time, and it has been found that the spectrum has a 
singularity of the Dyson type [S] in which spectral properties in the vicinity of the band 
center are functions of \n(E) where E is the energy relative to the center of the band. 
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The question of the band center localization has also been argued at length, and it 
was claimed that the band center mode may be weakly localized because of the strong 
fluctuations [TTlll2flU] . Ziman [T3"l[Ti] analyzing various uncorrelated distributions of 
the random coupling found the low-energy behavior of the integrated density of states 
(IDOS) 

1 V u 

IDOS-- = (1.1) 

2 \ln{E 2 )\ 2 V ' 



and the inverse localization length {ILL) 

ILL = W\< (L2) 



in which E is the energy and the logarithmic variance [T3"1[Ti] V u = | ((ln( J 2 )) 2 ) — (ln( J 2 )) 2 
is the amplitude factor, J being the off-diagonal coupling specified below. For defini- 
tions of the ILL and the I DOS see Ziman's work [13J particularly his discussion after 
the equation of the motion. Notice that the averages above are independent of the 
lattice site n since uncorrelated random variables are involved; otherwise they are site 
dependent. Concerning the shape of the IDOS and the ILL, according to Ziman's 
analysis [13J, they should depend on disorder only through a multiplicative constant, 
the variance in the distribution of ln(J 2 ). As a result, the spectra near E = have 
the same form for all sets of independent, identically distributed random variables 
J. In this context, such a property is referred to as scaling. When there is scaling, 
the data can be collapsed onto a single curve by dividing through by a factor that 
depends on the properties of the disorder but is independent of the energy. Although 
the above mentioned results were developed for uncorrelated disorder, our analysis 
shows that qualitatively similar scaling behavior for limited range of concentrations 
and energy can be obtained for a system with exponential correlations but with a 
modified amplitude factor V . Moreover, for a certain range of correlation lengths 
a higher order expansion of the Dyson singularity can reproduce the numerical data 
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better than just a single term. 

A disordered ID exciton model with nearest-neighbor interactions has the following 
equation of motion 

Jn+\U n +l + Jn-lUn-1 = EU n , (1-3) 

where U n denotes the usual exciton creation operator and J n are the couplings. This is 
equivalent to the XY model in a strong applied field in the z-direction (< Sz >= S). In 
our numerics we calculated iteratively the Lyapunov exponent [T3][T4] 'y(E) — In jf- 
(see appendix), where N is the number of sites [T3p4] . whose real part is related to 
the ILL while the imaginary is related to the IDOS. We consider couplings that 
have correlated or uncorrelated (for comparison) disorder and are of a dipolar type. 
Since the displacements of the atoms in quasi-one dimensional systems are likely to 
be small compared to the lattice constants, it is usually enough to consider the first 
order term in the expansion of J in terms of the displacements [15]. For the first 
order expanded coupling we have J n = 1 + A^ n with unit lattice constant where A 
is the strength of the disorder and £ n is a random sign variable that is exponentially 

| m — n\ 

correlated < £ n £ m >= e '( c ) where 1(c) is the correlation length as a function of 
"wrong sign" concentration c (see below). Thus the dipolar coupling J n has a non 
fluctuating part and a fluctuating part as in Dirac fermions with random-varying 
mass 08]. 

The generation of random binary sequences with correlation was deviced by imple- 
menting the convolution method of constructing a random squence. Such squence can 
be used in various applications as in designing low-dimensional devices with desired 
properties [7ffl6] . A detailed discussion on this topic recently appeared [17]. How- 
ever, a more efficient and easily handled correlated sequence, related to the spin 
glass chain problem [TECH], can be obtained through uncorrelated random num- 
bers as outlined below. The correlated random number at site n follows the relation 
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£ n = t^n-iXn = X1X2X3 . . . x n -\x ri with the distribution (of Xi i within [1, n]) 



P{xi) = (1 - c)5(x t - 1) + c5(xi + 1), 



(1.4) 



The Xi are uncorrected between different sites and c is the concentration of "wrong 
signs" such that c = and c = 1 mean no disorder. Clearly £ n is exponentially 



— l/ln\l — 2c|. For the numerical computation uncorrelated Xj are obtained from 
the uniformly distributed random number generator sequence. For the uncorrelated 
results we set £ n = x n at the particular concentration c = 0.5 where both correlated 
and uncorrelated cases produce identical results since Z(0.5) = 0. 

It should be emphasized that while the model for uncorrelated disorder and the model 
for correlated disorder involve the parameter c, in the latter model c determines both 
the fraction of wrong sign bonds and the correlation length. The two models intersect 
at c = 0.5 where there are equal numbers of ± interactions and the correlation length 
is 0. We have carried out extensive numerical studies to verify the equation for the 
correlation length given above. 

We note here that unlike the off-diagonal case above, many works have recently ap- 
peared on the diagonal tight binding model with correlated disorder dealing with 
various aspects of the problem [2"D"|21f22| . It was revealed that even short range cor- 
related disorder at the band edge caused the anomalies in the spectrum [23J. Frenkel 
excitons in off-diagonal and diagonal models have markedly different spectra, i.e., 
the former has a Dysonian singularity [T3l[T4"] at the band center while the latter 
has a power-law singularity at the band edge. Hence the short-range correlations are 
expected to influence off-diagonal spectra more than in the diagonal case. 



correlated (£ n £m) = (1 — 2c) 



where one can define a correlation length [24| 1(c) 
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Fig. 1. The integrated density of states (IDOS, left panel) and the inverse localization length 
(ILL, right panel) plots for the correlated distributions of the dipolar type couplings with 
A = 0.25, 0.50, 0.75 bottom to top. Symbols are data for a chain of 10 8 sites for respective 
concentrations divided by f(c) = The lines, shown as guides to the eye, are plots of 
the amplitude function /(0.5)(\/2 tanh -1 A) 2 divided by ln(E 2 ) for the ILL or its square 
for the IDOS. The parameter c is related to the correlation length 1(c) by the equation 
1(c) = -l/m|l - 2c|. 

2 Results 

The distribution and localization of the exciton modes are characterized by the Lya- 
punov exponent. In the case of one-dimensional arrays with nearest-neighbor inter- 
actions, the Lyapunov exponent can be determined by making use of mode-counting 
techniques [25126] . Recently we studied [27] similar excitons with orientationally dis- 
ordered couplings in conjugated polymers and calculated the full spectrum and optical 
line shapes. The goal here is to look at the effects of exponentially correlated disorder 
on the behavior of the IDOS and ILL near the center of the band (E = 0). We have 
produced three sets of data for each case varying the correlation length i.e., c) and 
the strength of the disorder, A, with A = 0.25, 0.50, 0.75. 

Figure [T] displays the results of the ILL and IDOS for exponentially correlated distri- 
butions of the couplings. For fixed A, we divided the data produced for a particular 
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concentration c by a scaling factor f(c) = c/(l — c) (see next section and the ap- 
pendix). The data for selected concentrations are shown in Fig. [1] for the indicated 
values of A. We expect first that for a range of concentrations the data should collapse 
to a Dysonian singularity for fixed A. However, Fig. [1] reveals that as c increases, the 
ILL and the IDOS take values that differ increasingly from predictions based on 
Ziman's analysis [32] and a slow approach to the band center is observed. In this 
limit, the data clearly show qualitative scaling behavior but deviations with varying 
magnitude exist as c — > 1 and c — > where the correlation length is infinite. Notice 
also that the deviations are stronger for the IDOS than for the ILL. 

Clearly as energy decreases, the scaling gets better; however, the deviations become 
worse as one gets closer to c — > 0, 1. To see whether longer chains and/or lower energies 
improve the situation, we checked further and found that the longer chain did not 
improve the agreement appreciably, but as the energy is lowered, the results follow the 
scaling behavior more closely up to a certain limiting c. Similar disagreement between 
the data and the Dysonian (with a single term) singularity were encountered in the 
SUSY [7] off-diagonal tight binding case for continuous exponential correlations. The 
analytical investigations [8] revealed that the spectra for the exponential correlation 
can be expanded for the IDOS as 

IDOS - ~ = a x /\ ln(£ 2 )| 2 + a 2 /\ln(E 2 ) | 3 + a 3 /\ln(E 2 )\ 4 + . . . , (2.1) 

and for the inverse localization length {ILL) 

ILL = b x /\ln{E 2 )\ + b 2 /\ln(E 2 )\ 2 + b 3 /\ln(E 2 )\ 3 (2.2) 

where oi, a 2 , a 3 , b\, b 2 are given as of functions of the correlation length multiplied by 
the strength [7]l8] of the exponential correlations. The equations fll.ip - fll.2l) are only 
the first order term in the expansion and a\ and b± should correspond to correlated 
amplitude factor V c (c). Figure [2] clearly shows the improved results obtained using 
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Fig. 2. The integrated density of states (IDOS, left panel) and the inverse localization 
length (ILL, right panel) plots for the correlated distributions of the dipolar type couplings 
with A = 0.50. Symbols are data for a chain of 10 s sites for particular concentrations (where 
significant deviations occur) c = 0.01,0.05,0.1,0.9,0.95 from bottom the top and the lines 
are the scaling function /(c) (\/2 tanh -1 0.5) 2 divided by ln(E 2 ) for the ILL or its square for 
the IDOS plus the higher terms given in Eqs. (|2.1P "( [2T2|1 where /(c) = j^—. The parameter 
c is related to the correlation length 1(c) by the equation 1(c) = — 1/ln |1 — 2c|. 



the Eqs. fl2.1l) - fl2.2p for the selected concentrations where the worst deviations followed 
from the first order result: c = 0.01, 0.05, 0.1, 0.9, 0.95. Unfortunately, for c < 0.01 and 
c > 0.95, this type of expansion appears to fail; however, other forms of fitting with 
fractional powers can give better results such as tried in the SUSY case [7118] . We will 
not here discus significance of this fractional fitting since there doesn't seem to be 
any analytic theory that predicts it. 



3 Discussion 



To shed some light on the scaling behavior for the correlated case, we first compute 
the binary uncorrelated case. A simple calculation of the logarithmic variance, the 
amplitude of the spectrum, for the uncorrelated case given in Eqs. (Il.ll) - fll.2p yields 
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the result 

V u (c) =4c(l-c)(v / 2tanh"M) 2 . 



(3.1) 



The amplitude factor for the uncorrelated disorder clearly displays a scaling behavior 
in concentration c as well as in A. We can identify the scaling functions for the con- 
centration g(c) = 4c(l — c) and the strength of the deviation h(A) = (a/2 tanh -1 A) 2 . 
The top inset in Fig. [3] displays the uncorrelated data vs concentrations. The three 
lines are the calculated amplitudes V (c) = h(A)g(c) presented in Eq. (13. ip for 
A = 0.25, 0.50, 0.75. As expected, the theory for the uncorrelated case works well 

The amplitude factors in Fig. [3] (large panel) are found from the best fit to the 
Dysonian singularity. They are obtained matching the smallest value of the energy 
rather than the value obtained by fitting over a range of energies. In this way, we will 
be more likely to get a value characterizing the asymptotic region. The calculated 
amplitude scaling function is also shown in Fig. [3] for the correlated case V c (c) = 
c(l — c) _1 ^v / 2tanh _1 Aj for the three specified A values above. In the appendix 
we will argue why this formula holds. We notice immediately that the concentration 
dependence has a different form f(c) = j£- as compare to g(c) = 4c(l — c) of 
the uncorrelated case while the scaling on A is the same for both cases. h(A) = 
("s/2 tanh -1 Aj . This asymptotic fit can further be checked by plotting the a x and b\ in 
the expansion above since they correspond the coefficient of the Dysonian singularity. 
The bottom inset displays the a\ and h\ obtained by the best fit as compare to 
the theoretical line V c (c) (for A = 0.5). Both plots clearly support the asymptotic 
amplitude scaling. 

Although the results obtained by asymptotic fit are clearly reproduced by the theo- 
retical amplitude factor, there are striking differences between correlated and uncor- 
related cases. What is particularly interesting and not expected is that the results 
for correlated disorder is not symmetric about c = 0.5 like those of the uncorrelated 
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Fig. 3. The asymptotic amplitude factors for correlated disorder vs concentration c with 
A = 0.25, 0.50, 0.75 bottom to top. In the main plot, the lines show the correlated am- 
plitude V°(c) = /(c)(\/2tanh _1 A) 2 where f(c) = j^. The top inset is the same plot 
for uncorrelated case for the same A values. The lines are a plots of the uncorrelated 
amplitude V u (c) = g(c) (\/2 tanh -1 A) 2 where g(c) = 4c(l — c). The bottom inset shows 
/(^(v^tanh -1 A) , the expansion coefficient of the IDOS and the ILL a\ and b\ , re- 
spectively vs c for A = 0.5. The parameter c is related to the correlation length 1(c) by the 
equation 1(c) = — 1/ln |1 — 2c|. 

case. Also, from our scaling analysis the ILL coefficient becomes large as c — ► 1, but 
vanishes as c — > 0. This is surprising since the modes for c = 1 are extended just 
like the modes for c = 0. The divergence of the amplitude factor as c — > 1, is not 
inconsistent with the results for c = 1 since the energy interval over which scaling 
holds shrinks to zero in the same limit. This can be seen from Eqs. fll.lj) — (jl.2p by 
setting either the IDOS or the ILL to a constant value and solving for the cut-off 
energy as as a function of c. As shown in the appendix, the amplitude factors perform 
a correlated random walk in parameter space. For our case, the square of the length 
of the random walk is proportional to cN/(l — c) not iV alone as in the uncorrelated 
case. That is the source of the marked difference between the two cases. As the con- 
centration approaches 1, divergent behavior occurs. However, we note that we cannot 
take c arbitrarily close to 1; beyond a certain point, the random walk sequence can 
not be constructed [i~4"] . 
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4 Summary 



We have investigated the dynamics of ID Frenkel exciton systems with correlated 
off-diagonal disorder. We have used a negative eigenvalue counting technique [25] 
which provides a simple and physically transparent analysis of the IDOS and the 
ILL. We investigated the question of scaling when there was correlated off-diagonal 
disorder such as might occur, for example, when the off-diagonal interaction depended 
on the distance between two ions and thus would be affected by small displacements 
of the ions from their equilibrium positions |15j . Our numerical data indicate that 
in the correlated case the scaling behavior found by Ziman [T3] is followed only for 
a limited range of c and very small E . We compared our results with the scaling 
predictions of Ziman [13] and with the similar tight binding correlated electronic case 
[Tf5] . We computed the asymptotic amplitude factors as a function of concentration 
and found that they are equal to the uncorrelated variance at c = 0.5 multiplied 
by the factor -^- c (see appendix). This factor [T8I28] has played a very interesting 
scaling role in number of unrelated problems. We found that with the new variance 
for correlated distributions rather unexpected behavior was obtained as compare to 
that of the uncorrelated case. In particular, the asymptotic amplitude factor is not 
symmetric about c = 0.5 while the ILL diverges for concentrations approaching 1 
in contrast to what is observed in the uncorrelated case. This amplitude factor can 
explain only the first order expansion in a Dyson-type singularity as supported by 
the SUSY model [7118] . As shown in the appendix, if the random walk observation 
of Eggarter et al. is implemented, with correlation, the obtained amplitude factor 
can reproduce the data within a shrinking energy range as E — > as c is increased. 
Finally, we should emphasize that our work is a numerical study supplemented by the 
approximate theoretical analysis outlined in the Appendix. As mentioned previously, 
the surprising result is the discovery that near the center of the exciton band, a model 
with correlated disorder showed asymptotic behavior that is similar to the behavior 
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of systems without correlations, the only difference being in the amplitude factor 
appearing in the limiting expressions for the IDOS and ILL. Like those of essentially 
all numerical studies, our results are approximate. We hope that the availability of 
the numerical findings together with our analytical results will stimulate rigorous 
analyses of the model that may shed light on the origin of the similarity. 
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6 Appendix 

In this section we present the amplitude scaling related to the first order expansion 
coefficients a± and b\ in Eqs. fl2.ip - fl2.2l) . Similar scaling but in different cases has 
appeared before [29l30|31fl8H28] . However the off-diagonal problem here is more in- 
volved than those cases. First the Lyapunov exponent at the band center E = is 
calculated using Eq. (11.31) that can be rearranged [5|13lll4j 

11 1 in 

7(0) = - HUn/U,) = - £ \n(U 2n+1 /U 2n ) = - £ HJ 2n+1 /J 2n ) + -. (6.1) 

JV iv n=l ly n=l Z 

The bipartite nature (the chiral symmetry) of the Eq. (11. 3p is responsible for the 
logarithmic variance as rigorously shown in [32j and eventually the Dyson- type 
singularity [H] since 7(0) executes a random walk with a step A n = m(J 2n +i) — 
ln(J 2n ). For the uncorrelated case the averages yield < A n >= and < A„ 2 >= 
|[< (ln(J 2 )) 2 > -< ln(J 2 ) > 2 ]. The Lyapunov exponent p^l32] takes this form 7(0) = 
^V /< (E„=i A n ) 2 > + f . The imaginary part here IDOS = = 0.5 and the 
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real part is related to the ILL. To make our point, we can attack the problem 
with a different angle. The average of the summation < J2n=i A n >= for both 
correlated and uncorrelated cases; whereas, the average of the squared summation 
J2 n =i < An 2 > +2J2 m ^n < A„A m >. It further reduces to N < A 2 > since the 
second sum is zero for the independent random numbers but non zero for the corre- 
lated random numbers. For the correlated disorder the summation is not easy. The 
complication arises when correlation present since < ln(J n ) > is not equal to its en- 
semble average rather it dependents on n. But this complication can be avoided in a 
following way. Let us calculate the Taylor expansion of ln(J n ) = ln(l + A£ n ) which 
reads 

ln(l + Ai n ) = A£ n - {Ai n f/2 + (Ae„) 3 /3 - (A^f/A + (A£ n ) 5 /5 + .... (6.2) 

As £ n 2k = 1 but for odd powers £ n 2fc+1 = £ n . Arranging even and odd powers we get 
ln(l + Ai n ) = (A + A 3 /3 + A 5 /5 + A 7 /7 + . . .)£„ - (A 2 /2 + A 4 /4 + A 6 /Q + ...). 



The first sum (with odd powers of A) is just the expansion of tanh~ 1 A and let 
the second sum (with even powers) be D which will be eliminated when it is in- 
serted in A n = tanh -1 A(^2n+i — 62™) ■ The same expansion can be used to calcu- 
late the uncorrelated case as well; the step takes the form A n = tanh^ 1 A(x2 n +i — 
X2n) (recall that number of steps N/2). Using Eq. (11.41) . the variance is given by 
V u (c) =< J: En A n 2 >= 2 < A 2 >= 4c(l - ^(V^tanh^A) 2 so that the ILL be- 
comes 5ft 7 (0) = ^V u {c)/N = \Z2tanh~ 1 A^J 4c(l - c)/N. For the correlated case, we 
have V c (c) = % < (En An) 2 >= {tanh^A f < (£„(6n+i - £ 2n )) 2 >. Taking the 
square, we get three terms < (£ ^odd/even) 2 > and twice < T,£oddY,£even >, hence 
the summation results in where inverse of this factor is encountered in previous 

works [28|18II29||3~1~] . The correlated random walk results in V c (c) = -^ c {\[2tanh~ x A) 2 
and the ILL = 3*7(0) = ^Jv c (c)/N = y/^tanh^AJ-^^. In contrast, the IDOS 
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value at the band center is not affected by the correlations since the distribution in 
eigenvalues is symmetric about the mid-point. The calculation above holds at the 
band center; however Eggarter et al, [H] have shown that for a certain range of E 
close to the band center, the random walk behavior still held and the ILL and I DOS 
can be represented as in Eqs. fl2.ip - fl2.2l) where the amplitude factor is the same as 
calculated at the band center. 
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